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We establish an intriguing connection between geometry and thermodynamics in 
the critical g-state Potts model on two-dimensional lattices, using the g-state bond- 
correlated percolation model (QBCPM) representation. We find that the number of 
clusters < A'^c > of the QBCPM has an energy- like singularity for q ^ 1, which is 
reached and supported by exact results, numerical simulation, and scaling arguments. 
We also establish that the finite-size correction to the number of bonds, < >, has 
no constant term and explains the divergence of related quantities as g ^ 4, the 
multicritical point. Similar analyses are applicable to a variety of other systems. 
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Percolation and the g-state Potts model (QPM) P are related to many interesting 
problems in mathematics and science and are ideal models for studying critical phenomena 
||T|-p. In recent years, much attention has been paid to universal quantities at or near the 
percolation point, such as the critical existence probability Ep or crossing probability [0, 
finite-size scaling functions excess cluster numbers etc. In a recent Letter, Ziff et 
al. @ calculated the number of clusters per lattice site, n, in percolation on two-dimensional 
lattices with N lattice sites and periodic boundary conditions (PBC). They found that 
n = ric + b/N + . . ., where Uc is n in the limit N ^ oo and 6 is a positive universal constant 
that may be calculated using conformal field theory (CFT) 0. In this paper, we consider 
the g-state bond-correlated percolation model (QBCPM) on planar lattices G of sites 
and E bonds, which is equivalent to the QPM on G; the numbers of bonds and clusters 
of a subgraph G' of G are denoted by Nh{G') and Nc{G'), respectively. In the QBCPM, 
as in ordinary percolation, a natural focus on geometric properties such as cluster number 
arises. However, the system also has non-trivial thermodynamics, impelling an investigation 
of the connections between geometry and thermal behavior. In this work we concentrate on 
the critical Potts models for definiteness; however, our methods are much more generally 
applicable, as pointed out below. 

We address this question by investigating the universal behavior of finite-size corrections 
(FSC). We show by exact calculation that when g 7^ 1, the FSC for < A^^ > is linearly 
related to the FSC for < >, i.e. surprisingly, the number of clusters has an energy-like 
singularity. This is quite different from the case g = 1, which is equivalent to bond random 
percolation |jTO| studied by Ziff et al. Numerical simulation, scaling theory for the 



infinite system, and finite-size scaling arguments verify and illuminate this conclusion. The 
latter also implies that < Ni, > has no constant finite-size scaling term at criticality, which 
we verify explicitly for the Ising model on a square (sq) lattice. We also find that the FSC 
of < A'^c + gNi) > and its higher cumulants (where g is defined below, and (7 = 1/2 for a sq 
lattice) diverge as g — 4, which is attributable to the onset of logarithmic corrections at the 
multicritical point and is understandable from a renormalization group (RG) picture. 



Here we briefly review the connection between the QBCPM and the QPM P,[T]|. In the 



QPM, each site of the lattice G is occupied by a spin Si with spin components —s, —s + 
1, s — 1, and s, where 1 < i < N,2s + 1 = q, and q is an integer. The Hamiltonian of the 
QPM is given by 

- H/ksT = KY1 ^{Sii s,) + Bj2si, (1) 

<i,j> i 

Here the first summation is a sum over all nearest neighbors, S{si, Sj) = 1 or when Sj = Sj 
or Si 7^ Sj, respectively, K = J/ksT > is the normalized NN coupling constant, and 
B = h/ksT is the normalized external field with being the Boltzmann constant and T 
being the absolute temperature. 

Using the subgraph expansion of Eq. (|l]), Hu has shown that phase transitions of the 
QPM are percolation transitions of the QBCPM, in which a subgraph G' appears with the 
weight 



TT 



where p = 1 — exp {—K) = 1; the spontaneous magnetization M and the magnetic suscepti- 
bility X of the QPM are related to the percolation probability P and the mean cluster size S 
of the QBCPM, respectively. These connections ensure that phase transitions of the QPM 
are percolation transitions of the QBCPM [Q. The partition function of the QPM at zero 
magnetic field may be written as 

= $:(exp (K) - lf^(^')q^^(G') = exp (KE) ^ 7r(G',p, q). (3) 

G' G' 



Here the sum is over all G' of G |]I0|. The internal energy U and the specific heat G^ of the 
QPM are related to the average number of occupied bond, p, and the fluctuations of the 
number of occupied bonds, G2b, of the QBCPM, respectively 0. 

Using the Swendsen-Wang algorithm |T^, we calculate the average number of clusters 
per site n of the critical QBCPM on L' x L square lattices with PBC in both horizontal 
and vertical directions; the number of spin components q is an input parameter taken to be 



1, 2, 3, and 4. It should be noted that n in the hmit L', L — > oo, denoted by nc, follows 



from exact results for the critical Potts free energy on several planar lattices P,p^. We 
plot — ric as a function of in Fig. 1 which shows that the data for q = 1 are on a 
linear curve. The linear least-square fit of these data gives rtc = 0.09807(6) and the slope 
b = 0.884 ± 0.002, which are consistent with the result of Ziff et al. However, results 

for q=2, 3, and 4 are quite different, namely the curves for q >2 have negative slopes, which 
suggests that the argument of Ziff et al. p to relate the slope b to the average number of 
clusters wrapping around the toroidal system is invalid and signals a new behavior as we 
show below. 

To understand the curves in Fig. 1 for q > 2, consider the partition function Z^. of the 
planar lattice QPM at the critical point Pc = 1 — e'^": 

Zc = umf'^'''Y^^'''^- (4) 

G' 

Here f{q) = e^^" — 1 and is known exactly for square, planar triangular, and honeycomb 
lattices 0; for sq lattice, /(g) = ^/g. Z^. is supposed to factor as Zc = ZnZu, where Z„ is a 
nonuniversal factor and the universal factor Zu gives FSC Exact results for Zu follow from 
the Coulomb gas formulas of Di Francesco, Saleur and Zuber (DFSZ) The cumulants 
Cn of Nc + gNb are given by C„ = [q{d/dq)]'^ In Zc, where g = g{q) = qf'{q)/f{q) and is 1/2 
for the sq lattice. Since Z^ = Zu{L'/L), FSC's to C„ are scale invariant. Thus, as in for 
q=l 

Cn = a^LL' + bn{L'/L) + 0(1/ L), (5) 



where 6„ is the universal FSC and may be derived from DFSZ [14]. It follows that there is 



no divergent FSC term for any C„ for g < 4. In particular, for n = 1 we have 

Ci =< iVe + gNb >= aiLL' + b{V / L) + 0(1/L). (6) 

For g = 2, 6(1) = 0.967734... and b{2) = 1.06463...; for g = 3, 6(1) = 1.05779... and 
6(2) = 1.13321 . . .. Since < A^"?, > is proportional to the internal energy, which has a 



singular FSC proportional to V-/'^ at criticality, Eq.(P) implies that the FSC for < > 
has an energy-like singularity with amplitude -g times the amplitude of < A'^^ >. A similar 
argument holds for any Cn, suggesting that Nc ~ —gNb in the sense of FSC, i.e. we can 
replace Nc by —gNb to calculate any leading FSC. 

This conclusion also follows from scaling for the infinite system. The singular part of 
the free energy per site fs may be written as fs ~ A{q)[p — Pc{(l)]'^~"^'^\ where A{1) = for 
(random) percolation and a is the specific heat exponent. Differentiating fg with respect to 
q, we find < N^ >~ —A.{q)p'c{q)[p — Pc{q)Y~^^'^^ LL' , showing that < Nc> is energy-like to 
leading order for q ^ 1. 

The universal (singular) part of the free energy is defined above at the critical point. 
According to finite-size scaling theory Fu also extends to large but finite systems near 



criticality, with 

F^ = LL'U^,lj{{[3-l3c)L^'nLL' 

^B + {/3- ^c)CL^'^ l^cfDL^'- + .... (7) 

Here /3c, i^-, B, C, and D depend on q; B, CL^^'^, and DL^^" determine the universal 
(singular) terms in the free energy, internal energy, and specific heat at the critical point, 
respectively. Let x = e^"^ — 1, then Xc = f{q)- Using the total partition function Z^r and 
< Nc >= q{d / dq)lnZ N , < N^ >= x{d / dx)lnZ ^ , we find for (3 = (3c that 

< Nc{G') >^ UcLL' - -llM_CL'/^ + qB'iq), (8) 

< Nb{G') >^ UbLL' + -li^CL'/r (9) 

J yXc + J- J 

Note that exact results for rib are available and Ub = 1 for the sq lattice Potts model for any 
q 0. Therefore, 

Ci =< Nc + gNb >= a,LL' + qB'{q), (10) 

where ai = Uc + gub, which agrees with Eq.(^ with b = qB\q). Note that Eq.(^ also 
implies that there is no constant FSC to < A";, >. 



It follows from finite-size scaling theory [|T5| that 



C2b =< N^>-< Nb >'= n2bLL' + caL^/" + . . . . 



(11) 



The CFT result therefore suggests that 



<N^> - <N, >2= n2cLL' + 9^021^^" + . . 



(12) 



< N^Nb > - < Nb>< >= UcbLL' - gc2L'^^'" + .... 



(13) 



It follows from Eqs. (|), ([TT|) , (|T2D , and ([T3|) that 02 = "^.20 + + '^g^cb- Now we proceed 
to test the above predictions. 



In 1 16], the internal energy of the Ising model on a large L' x L square lattice at the 
critical point, Ui{Tc)/Jj, is given by -Uj{Tc)/Jj = + 20/1 + d' / L"^ + . . ., where 9 = 
^2^3^4/(^2 + 6*3 + 6*4), Ji is the coupling constant of Ising spins and is related to J by 



J = 2Jj, 62, 63, and 64 are elliptic theta functions defined by Eq.(3.14) of [16] and d' was 
not determined in [|1^. We have extended the expansion of —Ui(Tc)/ Ji up to order 
and find 



(14) 



where 



e 



926364 



and Bi 



'K^R6l + 6l + 6l 



(15) 



^2 + ^3 + ^4 96 ^2 + ^3 + ^4 

The aspect ratio (i?) and the modulus {k) of the complete elliptic integrals of the first 
kind {K{k)) are related to each other by i? = K(k')/K(k) with k' = \J\ — k'^. 

Equation ([T^ ) shows that d! is zero, as predicted by Eq.(|D. Since Uj(Tc) = 2 J/ + U and 



(16) 



U = -j^ IhZn/N = p, p =< Nb > /E and n =< Nc > /N are given by 



^=2 



1 pcQ , Pc 



2L ' 2L3-ee. + oii3 



n = rir — 



p_,e 

2 L 



(17) 



where < Nb > and < Nc > satisfy Eq.(|^). As another test of Eq.(|^), we plot our p data for 
square lattice three-state Potts model as a function of L^/'^"^ = L~°'^ in Fig. 2, which shows 
that the data fit a linear curve with slope s = 0.1273 ±0.0005. In Fig. 3, we plot n — ric data 
for the Ising model and three-state Potts model as a function of 1/L for R = L'/L = 1 and 
2. The solid lines represent Eq. (|T^. The dotted line represents n — ric = —s/LP'^ + 
with h = 1.05779... obtained via [Q. The agreement between numerical data and our 
predictions is very good. 

From and the connection between the specific heat and the bond fiuctuations, C2b, 
of the QBCPM 0, we find that at the critical point pc of the Ising model C2b = C2b/LL' = 
C2 In L + n2b + 0[^). Here C2 = 2pl/7c = 0.218453 . . ., n2b = iP?5(0, R)/Kj + 1/{V2 + 1), 



and B{0,R) is defined by Eq.(4.21) of p[; for R = 1, n2b = 0.475235 . . .. Let C2c = (< 
> - < >^)/LL' and Cbc = (< N^Nb > - < Nb >< >)/LL'. For the Ising model 
(three-state Potts model), we fit C2b, C2c, and Cbc as linear functions of InL (L^/'^^^ = L'^-^) 
to obtain n2b, n2c, and Ubc and slopes. C2b — n2b, C2c — n2c, and Cbc — ribc for the L x L Ising 
model as a function of InL are shown in Fig. 4(a). The numerical values of C2 and n2b are 
0.21(8) and 0.47(6), respectively, which are consistent with exact values. The slopes for C2c 
and Cbc are 0.06(0) and -0.11(5), respectively, which are consistent with Eqs.(0) and (P^. 
C2b — n2b, C2c — n2c, and Cbc — Ubc for three-state Potts model as a function of L^/"""^ = are 
shown in Fig. 4(b); the slopes of these curves are 0.64(3), 0.16(6), and -0.32(7), respectively, 
which are also consistent with Eqs.(|rT])-(|l3]). 

As g — 4, the system approaches a multicritical point. From an RG point of view, its 
singular behavior may be understood in terms of a dilution field ip and temperature field 
r?| . Since ~ e = (4 — g) 5 , it follows from scaling theory that F„ will have an expansion in 



terms with integer powers of e along the line of critical points. Thus b, which is proportional 
to the q derivative of F^, and all higher cumulants 6„ diverge as g — > 4. This agrees with 



the results of a direct calculation using |]T3[, including the correct e dependence. For the 
cylinder geometry, b is finite but bn diverges for n >2, which is attributable to the vanishing 
of the leading term in the expansion of F„ in this geometry. 



For q = 4:, one cannot derive results for the FSC to < Nc > by differentiation. However, 



extending the scahng calculation in |TH], we find that to leading order 



p=<Nb> /E^ 0.5 + Ax{l - 2alna;)~^/^ = 0.5 + w{x), (18) 

where x = L^/*^^^ with z/ = 2/3 for four-state Potts model 0, A and a are non- universal 
constants and the sq lattice bulk value p = 0.5 has been used. The FSC part of this result 



includes the effects of the constant term in the scaling relation for the free energy [|T8[. In 



Fig. 5, we plot data of p — 0.5 and n — ric for four-state Potts model ||T3| as a function of 
X = Fitting p- 0.5 to w{x) of Eq.® gives A = 0.17±0.01 and a = 0.41 ±0.05. The 

solid and dotted lines in Fig. 5 represent w{x) and —w{x), respectively. Since E = 2LL' 
on the sq lattice, Fig. 5 shows that —w{x) also gives the leading FSC to n — ric, which is 
similar to the cases q = 2 and 3, i.e. we have numerical evidence for the relation Nc ~ —gNb 
when g = 4. 

Besides the Potts model, cluster representations are also useful for understanding critical 
properties of a model of hydrogen bonding in water, a dilute Potts model, the 0{n) model, 
quantum spin models, and many others P,p!^. Our methods are useful for understanding 
finite-size corrections in these systems. 

We are indebted to I. Affeck, A. Aharony, J. L. Cardy, M. E. Fisher, E. V. Ivashkevich, I. 
Peschel, P. Upton, F. Y. Wu and R. M. Ziff for useful discussions. This work was supported 
in part by the National Science Council of the Republic of China (Taiwan) under grant 
number NSC 88-21 12-M-OOl-Oll. 
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FIGURES 

FIG. 1. n — nc as a function of 1/Lp' for the QBCPM on L x L sq lattices with pbc (torus) for 
g = 1, 2, 3, and 4. 

FIG. 2. Numerical p of sq lattice three-state Potts model as a function of L^l^""^ with = 5/6 
for three-state Potts model. The solid line represents p = 0.5 + s/ L^'^ with s = 0.127(3). 

FIG. 3. Numerical n — nc of the L' x L sq lattice Ising model and three-state Potts model as 
a function of 1/L. The solid line represents Eq.(p!7|) for the Ising model with 6(1) = 0.967734. . . 
and 6(2) = 1.06463. . .. The dotted line represent the equation n — ric = —s/L^'^ + 6/L^ for the 
three-state Potts model with s = 0.127(3) and 6 = 1.05779 .... 

FIG. 4. (a) C2b — n2b, c^c — ?^2c, and c^c — ?^bc for the Ising model as a function of InL^ (b) 
C2fe — ^^26) C2c — ?^2c) and c^c — n^c for the three-state Potts model as a function of l?!^^"^ = L^-^. 

FIG. 5. p— 0.5 and n — ric for the four-state Potts model as a function of x = The solid 

and dotted lines represents w^x) and — w(x), respectively. 
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Fig. 1. Hu et al. 




Fig. 2 Hu et al. 
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